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We evaluate the modifications to the CMB anisotropy spectrum that result from a semiclassical 
expansion of the Wheeler-DeWitt equation. Recently, such an investigation in the case of a real 
scalar field coupled to gravity, has led to the prediction that the power at large scales is suppressed. 
We make here a more general analysis and show that there is an ambiguity in the choice of solution 
to the equations describing the quantum gravitational effects. Whereas one of the two solutions 
describes a suppression of power, the other one describes an enhancement. We investigate possible 
criteria for an appropriate choice of solution. The absolute value of the correction term is in both 
cases of the same order and currently not observable. We also obtain detailed formulae for arbitrary 
values of a complex parameter occurring in the general solution of the nonlinear equations of the 
model. We finally discuss the modification of the spectral index connected with the power spectrum 
and comment on the possibility of a quantum-gravity induced unitarity violation. 

PACS numbers: 04.60.Ds 



I. INTRODUCTION 

It is well known that black body radiation played a role 
not only in the historical development of quantum theory, 
but also in the formation of the global picture of modern 
fundamental physics, since both the observed microwave 
background radiation (hereafter CMB) in the Universe 
[l|-[5[ and the radiation from black holes predicted on 
theoretical ground by Hawking Q have an (approximate) 
black body spectrum. The former, discovered by Penzias 
and Wilson |2|, has revolutionized our understanding of 
cosmology, and the investigation of its anisotropics has 
shed new light on the physics of the very early universe 
thanks to the findings of the satellite missions COBE Q 
and WMAP [3] and other projects @. 

The other branch of modern fundamental physics that 
has motivated our research is the attempt of building a 
quantum theory of gravity 0, H|- This was seen for a 
long time as a logical step, independent of the ability 



* E-mail: binid@i cra.itl 
^E-mail: gcsposit@na.infn.it 
-^E-mail: kicfcr@thp.uni-kocln.de 
§E-mail: mk@thp.uni-kocln.de 

* * E-mail" ires . pcssina@gm ail.com| 



of performing observations: since gravity couples to the 
energy-momentum tensor of matter, and matter fields 
have a quantum nature, the scheme where the other fields 
are quantized whereas gravity remains classical can only 
have approximate validity. Although we cannot yet test 
physics at the Planck scale, we expect it should involve 
a quantum version of gravitation, so that both geome- 
try and matter fields are quantized. This area is where 
quantum theory and gravitational physics meet, leading 
to the unification of guiding principles as well as funda- 
mental interactions, and it involves the length scale out 
of which the present universe evolved. 

For a long time, it was thought that quantum gravita- 
tional effects, even when computable in a very accurate 
way, can hardly be checked against observations. Over 
the last two decades, however, attempts were made to 
establish a phenomenology for quantum gravity (see e.g. 
@)- One particular approach focuses on quantum grav- 
itational corrections to the functional Schrodinger equa- 
tion flol - [Ii| , as they are found from the Wheeler-DeWitt 
equation of canonical quantum gravity Q ■ This will also 
be the subject of this paper. An alternative canonical 
version is loop quantum gravity. In loop quantum cos- 
mology, analytic formulae for the power spectra of scalar 
and tensor perturbations suitable for comparison with 
observations were obtained [H|,[l6j]; alternatively, one can 
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explore the pre-inflationary dynamics [171 ] . 

Our paper is organized as follows. A brief review of 
fluctuations in quantum cosmology is performed in Sec. 
II, arriving at a coupled set of nonlinear differential equa- 
tions. Section III presents the solution of this system 
with and without a mass term, while quantum gravita- 
tional corrections are considered in Sec. IV. Section V 
studies in detail the possible predictions of enhancement 
or suppression of power at large scales, and Sec. VI is 
devoted to the spectral index. Concluding remarks and 
open problems are discussed in Sec. VII, while the Ap- 
pendix studies in detail the issue of possible violations of 
unitarity and how to get rid of them. We use units with 
H = c = 1 and a redefined Planck mass mp = y/3iv/2G. 



II. FLUCTUATIONS IN QUANTUM 
COSMOLOGY 



The treatment of small quantum fluctuations on a 
quantum Fricdmann-Lcmaitrc-Robertson- Walker back- 
ground was presented by Halliwell and Hawking fl8( | . 
These authors derived effective Schrodingcr equations 
for the various modes, which can be treated indepen- 
dently as long as the fluctuations remain small. In 
canonical quantum gravity, this corresponds to a Born- 
Oppcnhcimer type of approximation [7J . The time in this 
Schrodingcr equation is a "Jcffrcys-Wentzcl-Kramers- 
Brillouin (JWKB) time" that is defined from the vari- 
ables of the Friedmann background (scale factor a and 
homogeneous field tfi). 

Using the general Born-Oppenheimer approach pre- 
sented in Ref . [Io| , the authors of Ref . [IH have applied 
this scheme to fluctuations in quantum cosmology and 
extended it to the next order in »ip 2 to cover the first 
quantum gravitational correction terms. It was applied 
to the anisotropy spectrum of the CMB in order to derive 
the modification caused by these terms. 

Let us summarize here the main features of the for- 
malism for the derivation of the Schrodinger equation. 
For more details, the reader is referred to [l3| and the 
references therein. We choose a massive scalar field 
4> coupled to gravity in a spatially flat Friedmann- 
Lemaitre-Robertson-Walkcr universe; its fluctuations 
are expanded into Fourier modes with wave vector k ac- 
cording to 



The full Wheeler-DeWitt equation reads as |l| 



(x,t) = ]TA.(i) c 



(2.1) 



(The scalar- field modes fk are often called 6<fik-) Strictly 
speaking, this equation should be replaced by an integral 
representation, but we assume here that quantization is 
performed in a large box, so that we can take the wave 
numbers to be discrete for simplicity. One should, how- 
ever, bear in mind that the discrete sums in this paper 
have to be replaced by integrals if the universe is infinite. 
The relation between the discrete and the continuous case 
is discussed, for example, in the Appendix of Ref. (l9j . 



n + ]T?4 



fc=i 



where the Hamiltonians T-Lk of the fluctuation modes are 
given by 



1 * 

n, - 



3P_ 

dfl 



+ (fc 2 c 4Q +m 2 c 6Q )/2 . 
For the solution, one makes the ansatz 

oo 

{/fc}£=i) =M<x,<t>)i[Ma,<f>,fk), (2.2) 
fe=i 

and performs the redefinition 

ipk(ai,(t>Jk) = ^o{a,(f>)^k{a,(f>Jk)- (2.3) 
On writing 

MaJk) = e iS ( a > f «\ (2.4) 

and performing the expansion 

S(a, f k ) = mpSo + Si+ m p 2 S 2 + ..., (2.5) 

one can derive equations at consecutive orders in mp. 
At zeroth order in the Planck mass, one writes the fc-th 
component ip^ of the wave function as the product of a 
prefactor 7(a) with an exponential containing the phase 
Si(a,fk) according to 

^°W fe )=7(a)e iSl(Q < /fc) . (2.6) 
The JWKB time parameter t is then defined by 



9_ = _ e - 3a dSod_^ 
dt da <9a' 



(2.7) 



and each 1/^ is found to obey a Schrodinger equation of 
the form 



For pure exponential inflation, one has t = aH. 
^From a Gaussian ansatz 

^\t,f k )=^ 0) (t)e-^ 0) ^, (2.9) 

one then finds for A/^° and a coupled set of nonlinear 
differential equations, 

!<>(t) = -Ac- 3 «<^<>), (2.10) 



(0) 



(2.8) 



Iflf (t) = ie" 3 
dt k w 

+k 2 c 4Ht + m 2 e eHt 



rtk 0) (t) 



(2.11) 
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The second equation has the form of a Riccati equation, 
which has wide applications in physics (2p| . 

The time parameter (|2.7j) is the standard Fricdmann 
time t appearing in the standard form of the Robertson- 
Walker line element. In spite of this particular choice, 
the whole formalism is covariant with respect to time 
reparametrizations. Using the general Hamiltonian for- 
malism for canonical gravity, one has in fact N~ l d/dt 
appearing in (|2.7I) instead of d/dt, where N is the lapse 
function (see, e.g. the detailed explanation of this fact 
in Sec. 5.4.2 of |7j). Consequently, we can use any time 
parameter we like (e.g. conformal time instead of Fried- 
mann time) and directly rewrite all the results in our 
paper in terms of the new time. 



tion 



III. 



SOLUTION OF THE NONLINEAR SYSTEM 
WITH AND WITHOUT MASS TERM 



In this section, we solve the system of equations (|2,10[) 
and (|2.1ip with and without mass term. We shall here 
be more general than in Ref. We introduce the 

variable 



Ha(t) 



(3.1) 



The form of (|2.1ip suggests defining the dimensionless 
quantity p := m/H. It is then possible to obtain the 
general solution in terms of one unknown parameter, here 
denoted by U\, and the Bessel functions J„ and Y v with 
order 



1 



v 



(3.2) 



This solution is similar in form to the solution of the 
Klein-Gordon equation for a massive scalar field in 
de Sitter space (see e.g. Ref. [2l[, or Eq. (51) of Ref. 
[22j] . or Sect. 8.3.2 in Ref. [23j]). The reason is the general 
connection between the solution of the classical field and 
the solution for the function appearing in the exponent 
of the Gaussian wave function. 

In scenarios of inflation, one assumes that m < H in or- 
der to get fluctuations at super-Hubble scales with quasi- 
constant amplitude H/2n. This yields a real value for v. 

In explicit form, the solution of (|2.11[) reads 



1 



ev^o + mo) 



iUiYu+tf) + — ((3C7-1 + 2U 1 v)Y v (£) 



-2£J I/+1 (0 + (3 + 2i/)J 1/ (0) 



(3.3) 



For the massless case p = 0, (|2.11[) has the general solu- 



n? J (0 



k 3 



(Ui cost; — sin£)i 



H 2 £[(C/ X + cos f + (17^-1) sin^ 
£.3 i I Ji — U] J 



H 2 



(3.4) 



This coincides with the massless limit of the solution 
(|3.3[) . because then v = 3/2 and standard formulae for 
the Bessel functions of half-odd order show the agreement 
of (O with (15a . 

In the following, we shall restrict attention to the mass- 
less case. The massive case can also be dealt with along 
the following lines, but is technically much more involved. 
In Ref. [l3[ , the boundary condition was chosen that the 
quantum state (|2.9p approaches the free Minkowski vac- 
uum for large wave numbers k — > 00. This means that 



for £ -> 00 one demands that fij^ (£) ~ k 3 /H 2 t; 
In the massive case, this state is known as the Bunch- 
Davies vacuum [2~H |. For this boundary condition, 
reduces to the solution presented in Eq. (9) of Ref. [R 
Formally, it is achieved by setting U\ equal to — i, and 
one has then 



T 2 C 2 



ka 2 . 



O 



(0) 



/v :! 



1 



(3.5) 



The general solution is certainly richer than this one, 
because it involves ratios of combinations of Bessel func- 
tions of £, rather than just ratios of polynomials in the £ 
variable. On writing U\ = for the general case, we 
can re-express the in Q3.4p in the form 



^(0 = ^7 



AB* 



HH\B\ 



where 



(3.6) 



(3.7) 



A = p + icr, B = 7 + iS, 
and we have defined 

p := c(cos(/3 + £) + cos(/3-£)) -2sin£ 

= 2(£cos/3cos£ — sin£), 
a := C(sin(/3 + 0+sin(/3-£)) 
= 2£sin/3cos£, 

7 := c[(cos(/3 + £)+cos(/3-0) 

+£(sin(/? + 0-sin(/3-0) 
+2(£cos£- sin£) 
= 2C, [cos/3(cos£ + £ sin£) — (sin£ — £ cos£)] , 

5 := c[(sin(/3 + £)+sin(/3-£)) 

-fCcos^ + O-cosGS-O) 
= 2Csin/3[cos£ + sin£]. (3. 
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FIG. 1: The various curves correspond to values of £ = 
0.5, 2, 5, 27r, 10 and are ordered so that on the left part of 
the figure (say f3 close to 0) £ increases. In other words, the 
lower curve corresponds to £ = 0.5, whereas the upper curve 
to C = 10. 



As in Ref. [13[, we address here the classical quantity 
(Tfc (t) that is related to the quantum mechanical variable 
fk(t) of (|2.1[) by the following expectation value taken in 
a Gaussian state |V>fc)i i-e. 



- (V-fclA 2 !^) 



2Reft fc (t)' 



(3.9) 



In particular, one then finds at the level of approximation 



connected with ijjj^' , 

4° ) (*) = ^(Kefif( t ))" i . 
From the transformation rule 

At ^ H di 



one then gets 



*2°(t) 



In our case, this yields 



(3.10) 



(3.11) 



(3.12) 



H 2 



_ 2V2n 2 H 2 
ki 



f( P S- 


7cr) y 


1 — 




7cr) \ 


U(7 2 - 






U(7 2 - 


^ 2 )y 



^C{C + 2tt cos f3) 



y/sin /3 VC 2 + 47r cos /? + 47T 2 



£=2tt 

(3.13) 



The last step follows because £(i e xit) = 27r at Hubblc- 
scale crossing 13j. Figure 1 shows a plot of the function 



whose absolute value is taken here: 

V?(C + 27tcos/3) 



/(C/3) 



Vsin /3 v/( 2 + 4tt cos /? + 4tt 2 ' 



(3.14) 



when ( € (0,10), /3 G (0,7r). This function exhibits 
no relative minima or maxima, and hence there are no 
'preferred' values of ( and (3 from a mathematical point of 
view. The result (|3 . 13[) can be used to obtain the power 
spectrum filj[l5| 



p(°'(Jt) 



where 12 1 



^fe(ienter) 

such that we get 



^2 |^fc(iontor)| 2 



4<7fc(i) 



3 </>(t) 



^ (0) (fc) = ^ 2 |/(C,/3)| 2 -^ 



(3.15) 



(3.16) 



(3.17) 



IV. QUANTUM GRAVITATIONAL 
CORRECTIONS 

Proceeding with the expansion (|2.5|) to the next order, 
which contains terms proportional to rrip 2 , we arrive at 
a quantum-gravity corrected Schrodinger equation of the 
form fij 



i-^ (1) 



(i) 



^3a 



2m^f 



tf\ (4-1) 



where 



e 6Q iJ 2 , 
1 



^ 4 



fc 2 e 4 " + m'e' 



2„6a 



(4.2) 
(4.3) 



As in Ref. [13|, we make the following Gaussian ansatz 
for the corrected wave functions ipk ■ 

tf ) =(< ) w+-P 2 ^ i) (t)) 



xexp 



jK°)(t) + mp 2 n«(t))/|w 



(4.4) 



Here, a second term describing a possible violation of 
unitarity has been neglected in (|4.1I) . Such a term was 
found in the general derivation [lOfand can be inter- 
preted as follows (see also Ref. |2Q|). The general 
Wheeler-DeWitt equation is of the Klein-Gordon form, 
not the Schrodinger form. Expanding the conserved 
Klein-Gordon type current in powers of Wp 2 , one arrives 
at order m P at the exact conservation of the Schrodinger 
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(i) 



current, and at order nip 2 at a violation of this conserva- this reads as follows (since e _Q = H£/k) 
tion by a term that corresponds to the unitarity- violating 
term neglected in (|4.4j) . 

One can estimate that in most situations the unitarity- 
violating term is negligible compared to the correction 
term in (|4.4j) (Toj |. But one can also adopt the following 
viewpoint. After an appropriate redefinition of the wave 



2i£ 



d4 



(7 2 + $ 2 ) 



(C + LD) 



function, unitarity can be achieved at order nip [27 1. 
Such a procedure was also applied in the context of the 
quantum-mechanical Klein-Gordon equation in an exter- 
nal gravitational field pjj . Independently of which argu- 
ment is used, we shall no longer consider the unitarity- 
violating term in the main body of our paper, but we 
refer the reader to the Appendix for detailed calculations 
aimed at further clarifying this crucial issue. 

Inserting the ansatz (|4.4[) into (J4.ll) . one arrives at 



+ 





3jf(gWg) 2 ) 
2Vm 2 p 



+o(A 4 ) 



which can be cast in the form 

2 



(4.5) 



(4.6) 



1=0 



with time-dependent coefficients A21. Setting, in partic- 
ular, the overall coefficient A 2 of f 2 to zero, one finds the 
first-order nonlinear equation [l3f 



x on*) 



(d°\t)) -W k (t) 



W(t) 

Eventually, on defining 

C := pS — 7(7, D := pj + aS, 



(4.7) 



e 



H (C 2 -D 2 + 2iCD) \ 



( 7 2 + £2)5 



(4.9) 



the desired solution fij^ is taken to vanish 



In Ref. [1 

at late times. This expresses the idea that quantum grav- 
itational corrections should tend to zero at large times, 
which is certainly in agreement with observations. In the 
next section, we shall discuss a subtlety that arises when 
solving (|4.9[) with this boundary condition. 

V. ENHANCEMENT OR SUPPRESSION OF 
POWER AT LARGE SCALES? 



Taking as in Ref. [13j the initial state to be the ground 
state of Minkowski spacctime, one has to make the choice 
C = 1 and = |7T. For this case, (|4.9p reads 



(i) 



2ig (i) , 3^ (2g -i) 



d£ (£ - i 



2* (£-i) 3 



(5.1) 



The corresponding solution of (|4.9p is then inserted into 
the formula that generalizes (|3. 12|) at the next-to-leading 
order Hal. 



*?>(*) 



H(, d 
71 



Re^ 0) + mp 2 Rcn[}> (Z) 



from which one gets 



|C fc | 



(5.2) 



(5.3) 



Eventually, the square of the coefficients Ck yields the 
corrected power spectrum (see below). 

On requiring fij^ (£ = 0) = 0, one finds from (|5.1[) the 



(4.8) following exact solution: 



nf'fc) = -I^ [1 ^ } 2 2 )e2] + ^' part )(o, (5.4) 

where 

(i, P art)^ _ 3 [-5 + 6(1 + ifl + 4Ei (1, 2(1 + ifl) e^+g - 4g 2 (l + tQ] 

I 

and Ei(a, z) denotes the exponential integral defined by Here, we have to address the special case 

Ei M :=r^d*. (5.6) Ei ( M),r(0,,)^ lW , (5.7) 
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FIG. 2: The quantity ln(Cfe) is plotted as a function of 
H 2 /(k a ml) (see Eq. (JEU). 



cf. Ref. [291 ] . Sec. 5.1. From the exact solution (|5.5[) , one 
hnds for the coefficients Ck defined in (|5.3[) the result 



1 - 



54.37 H 2 



k 3 



7.98 ff 2 

fc 3 Ulp 



(5.8) 



see Figure 2. They correspond to an enhancement of the 
power spectrum compared to its value when quantum- 
gravity effects are neglected, 



V (1) {k) = 
~7>(°>(fc) 



89.54 H 2 



1 



fc 3 



fc 6 



ml 



-i 2 



(5.9) 



In Ref. [13|, another exact solution of (|5 . 1 [) is cho- 
sen (although it is there not given in explicit form). The 
solution has the same form as in (|5.5p , but with the expo- 
nential integral replaced by the other exponential integral 



(5.10) 



Ei(z) = - / — dt. 



f 



;,From this solution, one finds the value Ck given in Ref. 

m, 



c k 



1 - 



43.56 H 2 



k :i 



1 - 



189.18 H 2 



fc 3 



(5.11) 



This solution for Ck approaches 1 at large k (i.e. for 
small scales), but decreases monotonically to for small 
k (i.e. for large scales). In contrast to (|5.8[) . this solution 
thus leads to a suppression of power at large scales, 



V (1) (k) = V (0) {k)Cl 
123.83 H 2 



r {0) (k) 



i 



fc 3 



k 6 



El 

4 

mi, 



2 



.(5.12) 



What is the difference between both solution, that is, 
between the different versions of the exponential integral? 
Both solutions Ei(z) and E\(z) assume the value zero if 



FIG. 3: Plot of the imaginary part of Ei(l, 2i£ + 2) = £i(2i£ + 
2. 



£ = 0. Whereas, however, the solution Ei(z) approaches 
this value continuously (see Figure 3) , the function Ei(z) 
makes a jump with size ir in its imaginary part (29|, in 
agreement with the property (29| 

E x (-x ± iO) = -Ei(x) T 

where as above E\{z) := ^^-di when |argz| < n. Im- 
posing continuity as a reasonable selection criterion for 
the solution, since our JWKB ansatz for the wave func- 
tion should be differcntiable, would entail the choice of 
E\(z) and would thus lead to the prediction of an en- 
hancement of power at large scales, unlike the prediction 
of suppression in Ref. [l3| . 

Let us now revert to Eq. (|5.1[) . Remarkably, by passing 
to the new variable 



z = 1 + i£, 
it can be written in the form 



(i) 



dz 



and the solution reads as 

b 2z 3 



(5.13) 



1 

I 3 
(5.14) 



Hz' 



8z 2 + lOz - 5 



z- 



+4e 2z Ei(l,2z)] 



(5.15) 



with Pi determined so that fl^ (1) = 0. The passage to 
a complex independent variable is therefore very conve- 
nient in the course of performing and checking our cal- 
culations. 

Such a solution can be studied graphically by intro- 
ducing the complex polar representation for z = pc l8 and 
substituting it into the definition of ol 1 (z). On defining 
the functions 



f p {6) :=Re ^\pc ie ) , g p (6) 



Im 



(5.16) 

one then finds the behavior displayed in Figures 4 and 5. 



7 



i I 

i i 
i i 
i i 
i t 

i \ 

\f\ 





V\ e // x 

\ \\ A// i 

\ \ /// i 



/ i 



Hi \ 

[I 

i ! 7T 7c 3 ji 5 jt 3 7t 7n 

1 1 T T ~T~ "T" "T" ~T" 

!! e 

H 
i! 



- — p=l — -p=2 p=3 p=4 



•p=l- - p=2 p=3 p=4 



FIG. 6: Plot of the real part of the solution of (|5.22[) when the 
amplitude p of the complex variable z = 1+ 1 equals 1, 2, 3, 4. 



FIG. 4: Plots of /i, / 2 , fa, U (i.e. fixed values of p = 1,2,3,4) 
as functions of 9. As p increases, the curves show two en- 
hanced peaks at around 9 = and 9 = 2n. 



l/M 
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FIG. 5: Plots of pi, p2, P3, P4 (i.e. fixed values of p = 1, 2, 3, 4) 
as functions of 9. As p increases the curves show two enhanced 
peaks at around 9 — and 9 = 27r. 



The reader might still be worried by the fact that the 
limit as t — > oo mentioned before and in Ref. (l3j is only 
of mathematical interest, because it would correspond 
to eternal inflation, which is physically inconsistent. To 
answer this question, we begin by noting, from (|3.1[) . that 



1^ 



dt, 



and hence 



C(T) 



-H(T-t ) 



(5-17) 



(5.18) 



If T denotes the duration of the inflationary stage, it is 
enough to impose the boundary condition 



where, from (|5.18l) . 



C(T) = eato)- 



(5.19) 



(5.20) 



For example, if e is taken to be 10~ 2 , 1CP 3 ... and e is 
set to 0, we still find full agreement with the numerical 
values in Eq. (|5.8|) . 

If one is instead interested in the deep quantum gravity 
regime, one has to consider values of T so small that 
T — t < in (|5.18| . It is then appropriate to consider 
the variable 



in terms of which Eq. (|5.14[) becomes 



-SI 



(i) 



2i 



-O 



(l) 



3 J% ~ 2) 

2^(1-^)2 



(5.21) 



(5.22) 



1 



|, one can 



On defining the complex variable z 

then plot Sir' as a function of such a z. To perform 
a comparison with Figures 4 and 5, we plot in Figures 
6 and 7 the real and imaginary part of the solution of 
Eq. (|5.22[) when the amplitude p takes the same values 
considered in Figures 4 and 5. The plots corresponding 
to p = 3, 4 are virtually indistinguishable. The result 
is not, by itself, enlightening, but it shows that even the 
deep quantum gravity regime can be further investigated, 
if necessary 

Let us consider finally the general equation (|4.9p . Since 
the dependence on Q and /3 complicates matters, we limit 
our considerations to its linearization around Q = 1 and 
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2 4 



' p=l p=2 p=3 ~ ' p=4 



FIG. 7: Plot of the imaginary part of the solution of (|5.22[) 
when the amplitude p of the complex variable z = 1 + 1 equals 



/3 = 3/2ir and to the special case [i = 0, that is, we look 
for solutions of the form 



^(0 = + (C - W(0 + 08 - 3/27r)Q^°(0 , 

(5.23) 

Substi- 
tution of (|5.23[) into (|4.9|) gives the following equations: 



,(1)6/ 



where fii stands for the solution found in [n 



d^ 1)a 


2i£ 


d£ 


(i-0 


d"i 1)h 


2i£ 


d£ 





j(i)a 1 i^(i sin ^ — 2 cos£ — 2icos 2 £sin£ + 2 cos 3 £) 
fc 4 (isin£ - 3cos£ - 4icos 2 £sin£ + cos 3 £)(* - £) 3 



^ 1 + — 4(i- 4 ■ ^ ' (5 - 24) 

where 

4(0 = 12£ 4 - 6£ 2 + 18i£ + 3 + 8c 2i «P! + 12Ei(l, 2(1 + iO)e 2(1+i5) , (5.25) 



and Pi = -3/8 - 3c 2 Ei(l, 2)/2. These equations can be 
studied numerically. Unfortunately, the task turns out to 
be technically very hard, but these equations are given 
here because their investigation might shed light on the 
consequences of choosing another vacuum for the very 
early universe. This is the subject of future investiga- 
tions. 



VI. OBSERVABILITY OF THE CORRECTIONS 

We have used the uncorrected Schrodinger equation 
(|2.8[) to arrive at expression (|3.13p . from which we can 
immediately obtain the power spectrum (|3.17D . 



p 4 



\<Kt)\ 



This corresponds - apart from a dimcnsionless constant, 
which is not relevant for the following discussion - to the 
standard power spectrum of scalar cosmological pertur- 
bations llal. 



G 



H 1 



(6.2) 



where we have used the first slow-roll parameter defined 
as 

€ ~ H 2 ~ H 2 



it. 



(6.3) 



As we have seen in (|5.9[) and (|5.12|) . the quantum- 
gravitationally corrected Schrodinger equation leads to a 
modification of the power spectrum by a correction func- 
tion Cfe, such that we can translate this modification also 
to the standard power spectrum in the following way: 



Along the lines of [15| . we write 



^WDw(^) 




where ^wdw(^) either takes the form 



r(k) 



179.09 / H 

Up 



(6.4) 



(6.5) 



(6.6) 



"WDWVV - u3 

which follows from (|5.9[) , or the form resulting from (|5.12[) 

247.68 ( H 



J WDW 



"i-ttI- 



(6.7) 
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The basic equations in the theory of the spectral index 
n„ and its running a s involve the slow-roll parameters 

m 



7? := 



H(j> 



and H 



There is thus a change of sign in front of <5 W dw 
respect to the discussion in (l5j | . that is, 



with 



and 



o., 



1 := 



d log k 



dlogP, 
d log k 



2-q - 4e 



°"WDW 



2(5e?7-4e 2 



" ) ' ao WDW ' 



(6.9) 



(6.10) 



where use has been made of the approximate formula [l5| 
d 1 d 



dlogfc H At 



(6.11) 



jointly with the equations of motion. 

Giving up our definition that k is dimensionless and 
reinserting a reference wavenumber, which can cither cor- 
respond to the largest observable scale, fc m i n ~ 1.4 x 
10 _4 Mpc~ 1 [HI], or to the pivot scale used in the 
WMAP9 analysis, k Q = 0.002 Mpc" 1 [U, HI, we now 
write fc/fc m i n or fc/fco, respectively, instead of k. Since the 
ratio H/mp has to be smaller than about 3.5 x 10~ 6 be- 
cause of the observational bound on the tensor-to-scalar 
ratio r < 0.11 for k = 0.002 Mpc' 1 from the Planck 2013 
results [3l], HH, we find that for k — > k/ko the absolute 
value of the quantum-gravitational correction is limited 
by 



|*wdw(*o)| <2-2x 10" 9 , I 

while with the replacement k 
thcr weakened: 



J WDW 



(k )\ < 3.0 x 10" 



k/k m - ln this limit is fur- 



|^wdw(M| 



< 7.5 x 10 



-13 



J WDW 



(Ml < 1-0 x 10 



-12 



The difference for 5^ DW compared to [15J stems from 
fact that our upper bound for the ratio H/mp is weaker 
than the upper bound derived in [l5| from a different 
assumption. Furthermore, we have used the exact values 
of (16. pp and (|6.7[) instead of the approximate value 10 3 
as in jl5j . 

By comparing the quantum-gravitational corrections 
to the spectral index n s and its running a s derived 
above with the values determined from the WMAP9 
data, n s = 0.9608 ± 0.0080 and a s = -0.023 ± 0.011 
(using the WMAPg+eCMB+BAO+ffo dataset in both 
cases) [13, and the 2013 results of the Planck mission, 
n s = 0.9603 ± 0.0073 and a s = -0.013 ± 0.009 (using 
additionally the WMAP polarization data in both cases) 
[32^ , we see that our corrections arc completely drowned 
out by the statistical uncertainty in the data. Further- 
more, one can already rule out that further improvements 



of the statistics of Planck data and future satellite mis- 
sions to measure the CMB anisotropics more precisely 
will push these corrections into observable regions, as 
the main source for statistical uncertainty on large scales 
in the anisotropy spectrum is cosmic variance, which is 
given in terms of spherical multipolcs I = 2fc/fc m ; n by 
(see e.g. [33[) 



Var Ps {£) 



2£+l 



5.12) 



Defining £q as the multipole corresponding to the pivot 
scale feo given above, £n « 29, the conclusion of the 
detailed discussion in [15j, which determines the region 
in the {£,T S {(-) /Vs^o)) plane that is affected by cosmic 
variance can essentially be carried over to the present 
discussion. The only difference is the sign change for 
^wdw(^) an d a very slight suppression of the correction 
if one does not approximate the pre-factors in (|6.6I) and 
(|6.7[) . The order of magnitude of the correction stays the 
same and therefore also the conclusion that the quantum- 
gravitational correction is entirely negligible compared to 
the error induced by cosmic variance remains. 

We can repeat, however, for &^r DW the analysis to de- 
termine an upper bound on the energy scale of inflation, 
which was presented in [l3j for 5 



WDW 



Instead of assum- 
ing that C\ has to be greater than about 0.95 in order to 
be compatible with the observation that the anisotropy 
spectrum deviates from a scale-invariant spectrum by less 
than 5 %, we introduce the upper bound that C\ has to 
be smaller than 1.05 at k/ko ~ 1. From equation (|6.5[) 
and (|6.6[) . we then immediately obtain the upper bound 



H < 1.67 x 10~ 2 top w 4.43 x 10 17 GeV. 



(6.13) 



For i5 WDW the upper bound is [T3[ 

H < 1.42 x 10~ 2 top w 3.76 x 10 17 GeV. (6.14) 

Both constraints are clearly weaker than the already ex- 
isting observational limits of about H < 10 15 GeV, but 
it reassures that the present approach is consistent with 
these limits. 

As a final remark we want to add that the general dis- 
cussion in II 51 ] on non-Gaussianities in the squeezed limit 



arising from quantum-gravitational corrections is not af- 
fected by our results in the present paper, and we there- 
fore have to conclude that, at the present state, one can- 
not expect to see any effect from Wheeler-De Witt quan- 
tum cosmology in the bispectrum. 



VII. CONCLUDING REMARKS 

We have studied the quantum gravitational correc- 
tions terms to the CMB anisotropy spectrum as they are 
found from a Born-Oppcnheimcr type of approximation 
from the Wheeler-DeWitt equation. We have, in partic- 
ular, discussed a subtlety concerning the central equa- 
tion (|5.ip . Imposing the boundary condition ftfjp — > 
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for £ — > (corresponding to the absence of quantum- 
gravity effects at late times), we have found that there 
exist two solutions. One of them, which is related to the 
exponential integral E\(z), is continuous in this limit, 
whereas the other one, which is related to the exponen- 
tial integral Ei(z) makes a jump in the imaginary part 
of size 7r. Both solutions are in accordance with the re- 
quirement that quantum-gravity effects are unobscrvablc 
at late times. The numerical corrections to the CMB 
spectrum are of the same order in both cases. There is, 
however, a qualitative difference. Whereas the continu- 
ous solutions leads to an enhancement of power at large 
scales, the discontinuous solution leads to a suppression. 
If we adopt continuity as a condition for the allowed so- 
lutions, bearing in mind that the JWKB ansatz for the 
wave function should be differentiable, we have to predict 
an enhancement of power, unlike [l3| where a suppression 
is predicted. 

So far, these corrections to the CMB power spectrum 
are too small to be observed. Nevertheless, their size is 
much bigger than corresponding corrections in laboratory 
situations. One can thus express the hope that it might 
eventually be possible to test them in a cosmological set- 
ting. 

At a field-theoretic level, an interesting issue is whether 
the various choices of amplitude £ and phase /3 in Sec. Ill 
and V can describe physically relevant choices of vacuum 
other than the Bunch-Davies vacuum. A further issue 
seems to be whether our quantum cosmological calcula- 
tions can be used to test the recent theoretical prediction 
of circles in the CMB [34|. These are subjects for future 
work. 

For a full-fledged investigation, one has to repeat our 
whole analysis by using the gauge-invariant formalism 
of cosmological perturbation theory [3]J HU . While this 
will be postponed to a later paper, we emphasize that the 
main qualitative features and open problems are clearly 
seen already in the more restricted analysis presented 
here. 

Modifications of the power spectrum at large scales 
can also emerge from other situations. Inflationary mod- 
els with f2mj£ 1 lead to an enhancement of power for low 
multipoles [37|, both for open and closed models. A sup- 
pression of power is predicted from the consideration of 
a self- interacting scalar field [Hj], a model of just-enough 
inflation 13911 , or the introduction of non-commutative ge- 
ometry [40( 1 . Yet other work 41) has studied the nonlo- 
cal entanglement of the Hubble volume with modes and 
domains beyond the horizon, finding that this induces 
a dipole and quadrupole contribution in the CMB. The 
great challenge is, of course, to distinguish these contri- 
butions in possible observations. 



and to Giovanni Venturi for bringing the work in Ref. |27| 
to his attention. C. K. thanks Dominik Schwarz for dis- 
cussions. M. K. acknowledges support from the Bonn- 
Cologne Graduate School of Physics and Astronomy. 



Appendix A: The unitarity violation issue 

We are here going to see how the unitarity issue is 
solved with an alternative approach that enables us to 
explicitly verify that no unitarity violations occur. For 
this purpose, relying upon Ref. (27| . we start from the 
equation (see also Sec. 5.4.1 in 0) 



f h 2 d 2 
y^da* 



mpV(a) + Hm *(a, 



= 0, 
(Al) 

where Hm denotes a Hamiltonian depending generi- 
cally on matter fields (f> whose mass is light, -C 
rrip. By virtue of this inequality, introducing the Born 
Oppcnhcimer (hereafter BO) approximation 



*(o,0) = ip(a) X (a,4>), 



(A2) 



we obtain 
1 

2m 2 , V ^ da 



■da 2 



+ ipVax + i>H M X = 0, 
(A3) 

with Vg = mpV(a). In our BO approximation or, in 
other terms, adiabatic approximation, we can neglect 
the term -§^X because it varies slowly with respect to 
the scale factor a. Hereafter, in order to agree with the 
notation of [42[, we will use the Dirac-like notation so 
that x{ a i4>) = ( a >^|x)- Thus, evaluation of (x| on the 
equation obtained from (|A3|) upon replacing x(a, 4>) with 
{a, (j>\x) yields 



1 d 2 
2mp da 2 



i> + V G ilJ + {H M )i) + \ ( X 



d_ 

i)a 



X 



^=0, 
da 



(A4) 

where we have assumed (x\x) = / d</>x*X = 1- Now we 
obtain Vctp from the latter, and substituting it into (|A3[) 
gives 



H M - (H 



M 



\x) 



1 dlntp 



Wp da 
d_ 
da 



d_ 

da 



X 



x)\x) 



0. (A5) 



We can obtain the same result by keeping all terms and 
doing the approximation after the calculation. In fact, 
from (|A3[) . if we evaluate M on the left-hand side and 
add and subtract tp x ((xl^lx)) 2 ! we get 
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2mp 



D 2 + V G 



(H 



Mi 



1> = ~ 



2mp 



(D 2 )^, (A6) 



where we have introduced the operators 

— +iA, D:= — - 
da da 



D 



LA, 



(A7) 



11 



A~-i( X 



d_ 

da 



(A8) 



Following [43j . the condition (xlj^lx) = can be seen 
as a 'gauge' choice and it implies no restriction by virtue 
of the 'gauge invariance' of the full wave function. We 
choose (xlgjjlx) = as Eq. (|A4[) suggests because the 
equation we require for the semiclassical theory of gravity 
is that in which classical gravity is driven by the mean 
energy of matter. 

Evaluation of (\\ on Eq. (| A6|) and insertion into Eq. 
(|ATj) leads to 



i){H M - (H M ))x + ^D^D X 



2mp 



4>{D 2 - (D 2 )) X - 
(A9) 

This equation becomes equal to (|A5[) if we define the new 
wave functions 



exp 



X := exp 
* := ^x 



Ada' 
Ada 1 



0, 



X, 



(A10) 

(All) 
(A12) 



At this stage, the BO approximation is imposed by ne- 
glecting the terms (D 2 ) and (D 2 — (D 2 )) that are related 
to fluctuations. 

The next step is to require that ip should be the JWKB 
solution of Eq. ([A"4]l 



1 

N 



exp 



(A13) 



where S is the classical action defined by the trajectories 
that solve Eq. (|A4[) . and 



N = 



dS 



da 



= \/mpa, 



(A14) 



as one can obtain from the JWKB method. Now, sub- 
stituting this relation in (|A5|) , we obtain an equation for 

Ix) 



H M - (H 



Ml 



10- -\--\x) 

mp da da 



0. 



(A15) 



The latter might be regarded as a Schrddinger equation 
for the matter wave function 

(h m - iti^j \ X s) = 0, (A16) 

where we have introduced the time derivative through 

(A17) 



m— -[——— 

dr} rap da da ' 



and we have defined the new wave function 

rv 



\Xs) ■= exp 



(H M )dT]' 



Ix)- 



(A18) 



So we note that at the zeroth stage of the JWKB ap- 
proximation one obtains the usual evolution equation for 
matter (Schrddinger equation or, in the field-theoretic 
case, Tomonaga-Schwinger equation) Q. 

One may now search for a possible violation of uni- 
tary evolution for this matter wave function. In order to 
check this, we perform all the approximations after the 
calculation. On considering all terms, (|A5[) becomes 



*L r 

2 



■;/) 



dr) 
dlogN 
da 



Xs) = 

D + l(D 2 -(D 2 )) 



|X S )-(A19) 



Now we can investigate a possible violation of unitarity 
from the relation 



(Xl (h m - 



{xs\§-xs)-c.c. 



h 2 dlogN 
m?, da 



D 



2mp 



(D 2 -(D 2 ))\ X )-C.C. 



0, 



(A20) 



because we have to neglect (D 2 ) and (D 2 — (D 2 )) in the 
BO approximation. Thus, unless one considers a non- 
Hermitian Hamiltonian operator, there is no violation of 
unitarity. 

Now we want to find a relation between our approach 
and this one. It is indeed possible to rewrite the BO 
approximation in the form 



= exp 



(A21) 



and expanding S in powers of mp 
nipSo + Si 

nip 
1 

Si (a, 



mlS + Si + \s 2 + O ( \ 
nip \ m p 

1 



(A22) 



m p S (a) H 5-02(0) 

mi* 



where we have decomposed $2(0, <j)] 
Thus, the wave functional becomes 



-772(0, 



02(0) + 772(0, < 



1 



exp 



A^(o) 
Nk(o,) exp 

tpxXx- 



1 



—S 

h hrnp 



(T 2 



iSi 



V]2 

hmp 



(A23) 



In order to obtain the desired equations, we have to sub- 
stitute these expressions for ipjc and xk i n (|A6|) and 
[9|) . From (| A6|) we have three equations proportional 
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to O(mp), 0(m p ), C(m p 2 ), respectively, 



\s' ( ? + V = 0, 



(A24) 



<^> -i»%^+i»f =0, (A25) 



2 [ 2 N 2 



- 2- 



fi 2 TV 



A" 



<if M >_ 2 = -yCD 2 )o, 



(A26) 



where a prime denotes a derivative with respect to a. 
Analogously, from (|A9|) one obtains 



Fg, - (H M ) Q + ms, 




to C(m P ) and C(m P 2 ) respectively (where we have de- 
fined the c-number H M by H m \xk) = HmIxk"))- By 
comparing terms C(m P ), C(mp 2 ) we have 



(A29) 



Q/f Q! 2 

1 1 



-S' r& + H M A - S' a' 2 + \h-f - -f- = 0, (A30) 
respectively, where we have expanded Hm = + + 



+ As one can see, (|2.8[) and (|4.1j) are identical 

to (|A29P and (|A30P if we choose H M = H k . To perform 
the calculation of the possible violation of unitarity, it is 



extremely useful to note that the condition (x\D\\) = 
becomes 



(X\D\X) = (X\^\X)=0, 



(XKs\-^\XKs) = 



N K , S S{, S ^ 



1 JV^- ft mpfi, ' 



0. 



(A31) 



In the same way, as in (|AT9|) , for the C(mp 2 ) we obtain 




(A32) 



ST 

h 2 " hn k 

Performing the same calculation as in (|A"20]l . we finally 
obtain the same result 



ifr-7^(XKs\XKs 



i.S', 



Xk-^-Xk + C.C. 
oa 



0. (A33) 



To summarize, an appropriate redefinition of the wave 
functions leads to a description without unitarity viola- 
tion. This is why we can safely neglect the corresponding 
term in the corrected Schrodinger equation (|4.1j) . It is, 
however, an open issue which of the wave functions (the 
original or the redefined one) is the relevant one in the 
sense of a quantum measurement process. 
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